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ABSTRACT 


The prebuckling state and the stability behaviour of thin-walled 
isotropic and orthotropic shells under axisymmetric loads (axial load 
and internal pressure) is investigated, applying a nonlinear shell theory 
(finite deformations). 

At first the unbuckled equilibrium state which is axisymmetric 
under the given load is described. The derived system of ordinary non- 
linear differential equations is solved by a difference method with 
approximation of higher order. For isotropic conical shells, an exact 
solution for the linearized system can be obtained (linear bending 
theory) . 

In the second part of this paper the buckling loads are calculated. 
The buckling (bifurcation) state is governed by an eigen-value problem, 
the eigen-values of which are the buckling loads and the eigen-functions 
of which are the additional stresses and deformations. The eigen- 
value problem is also solved iteratively by the difference method. 

The influence of shell parameters and of an axisymmetric predefor- 
mation on the buckling load is discussed. The results are in good 
agreement with the experiments, particularly for a conical shell under 
axial load and high internal pressure. 
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1. 1 Introduction 


The stressed structures in aviation and rocket construction 
consist to a large part of thin-wall isotropic, orthotropic, and rein* 
forced shells. Detailed knowledge of the carrying capacity of thin- 
wall shells is required if the construction elements are to be designed 
as light as possible in favor of a high payload. 

The elastic stability behavior of cylinder shells has been intensively 
investigated during the last years according to the linear as well as to 
the nonlinear theory of stability. For conical shells only buckling loads 
have been determined so far according to the linear theory, because the 
basic equations describing the stress and deformation state of the coni- 
cal shell are more complicated compared with those for the cylinder 
shell. These loads have been determined among others by Mushtari 
and Sachenkov^, Schnell^ , and Seide^^»^^>^^ for the isotropic conical 
shell and by Leyko^^ , Serpico^^ , and Singer^'^^’^^ for the orthotropic 
shell. 

By means of the linear theory the qualitative effect of the geometri- 
cal characteristics and of the rigidity parameter on the buckling load 
can be determined. However, the resulting buckling loads deviate 
greatly from the buckling loads obtained experimentally, especially in 
the case, herein considered, of the conical shell under axial load and 
inner pressure. This discrepancy between experiment and theory is 
greatly due to the fact that according to the linear theory of stability 
the basic state of the stresses of the shell--that is the state of the shell 
before reaching the critical load--is replaced with the state of diaphragm 
stress, which neither satisfies the compatibility condition nor all 
boundary conditions (no circumferential expansion). 

In the following discussion, first the basic state of the stresses of 
the orthotropic conical shell (with the special case *'isotropic”) shall 
be described, taking into account nonlinear members. The derived 
nonlinear differential equations are approximately solved with the dif- 
ference method. With the given boundary conditions and for the load 
investigated (axial load and inner pressure), the solution is axially 
symmetric. 

For the linearized basic equations of the isotropic cone describing 
the state of the stresses and of the deformation under small loads at 
the edge, a general solution in closed form can be stated. 


1 


The second part of the discussion deals with the determination of 
the buckling loads. After exceeding the critical load, the axially 
symmetric basic state of the stresses becomes unstable and the shell 
buckles into a nonaxially symmetric state of equilibrium. The deter- 
mination of the additional stresses and deformations occurring in the f 

instance of buckling leads to an eigen- value problem, the eigen-value 
of which is the critical load. The iterative solution of this eigen-value 
problem is obtained by means of the difference method. 

If the calculation of the buckling loads is based on the basic state 
of the stresses which results from the basic nonlinear equations when 
the boundary conditions are completely satisfied, the buckling loads 
determined are in good agreement with the experimental values for high 
inner pressure. The discrepancy between the buckling loads calculated 
for an ideal cone and the experimental buckling values for low inner 
pressure (especially for p = 0) is essentially traced to the existence of 
prebuckling in the test shells. 

In order to be able to investigate the effect of such a predeforma- 
tion on the stability, the differential equations of the prebuckled conical 
shell are derived and iteratively solved. The applied method of solving 
permits us to take into account any arbitrary axially symmetric pre- 
buckling. 

The investigations described here principally apply also to the 
conical shell under outer pressure. 

1. 2 Glossary of Symbols 


bA 

B{y, <Po) 

^ik* ^ik 

Dik 

E 

Fj , Fin 

G 

h 

^k 

L(y) 

m 


Width of prebuckling 
Differential operator 

Stretching and shearing rigidity parameter 
Flexural stiffness and torsional stiffness 
Modulus of elasticity 
Stress functions 
Modulus of rigidity 

Distance between the equidistant intermediate points 
Height of the cone 
Differential operator 

Niimber of the equidistant intermediate points 
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m 




Mg, 

Me, 

Mge 

Ms, 

Me, 

Mse 

n 




Ns, 

Ne. 

Ns 

50 

Ns, 

Ne, 

Ns 

i0 

No 




Nv 




qo , 

Qo 




Qs> Qq 

P 

P 

P 

ri, T2 
s 

s 1 , S2 

t 

u 

V 

w 

WA 

z 


Number of bucklings. along the generatrix 
Section momentum 
Edge momentiam 

Number of bucklings in the direction of circumference 
Section forces 
Edge forces 

^2 2 
Si e 

Section force parcimeter No = — N„q 

Dn sw 

-No 

Dimensionless loadparameter N^^ = — — 
Deformation parameter of the basic state of stresses 



Transversal edge forces 
Inner pressure 

Dimensionless inner pressure parameter 



Axial load 

End radii of the cone 

Coordinate along the generatrix 

Length of the generatrix to the ends of the cone 

Wall thickness of the shell 

Local displacement of a point along the generatrix 

Local displacement of a point in the direction of 
the circumference 

Local displacement of a point perpendicular with 
respect to the central surface of the shell 

Predeformation perpendicular with respect to the 
central surface of the shell 

Dimensionless coordinate along the generatrix 
^ s 

z = £n-— 

Si 
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Zo 


= Hn 


S2 


Si 


ZA 


a 


Position of the prebuckling maximum along the 
generatrix 

Half-angle of the cone 


\s 0 
® s > ® 0 

0 , 0 

k 

^0> 

X 

V 


Shearing of a shell element s d 0 ds 

Stretchings of the central surface of the shell 
n 

sin a 

Circumference coordinates 0 = 0 sin a 
mk 
zo 


Curvatures of the central surface of the shell 

Sj cot a 


Shell parameter X = 
Pois son* s ratio 


T 


11 ^2Z 


n 


Potential energy of the deformed shell 


Index ( )o for stresses and deformations designates the state of the 
basic stresses 


Index ( )i designates the additional stresses and deformations, 
occurring in the instance of buckling. 


2. Basic Equations of the Orthotropic Conical Shell 


2. 1 The Law of Elasticity of the Orthotropic Shell 


The geometrical designations used in the investigations are 
given in Figure 1, Each point of the central surface of the shell can 
be fixed by the distance s and the circumferential angle 0 ; s and the 
angle 0 --where 0 = 0 sin a--represent the polar coordinate system of 
the development of the shell. The arrows of the coordinate axes point 
in the direction of positive values. Because of the singularities 
occurring at the tip of an axially loaded cone only truncated cones are 
investigated. The geometrical configuration of a truncated cone is 
uniquely determined by a, Si, S 2 . 

The law of elasticity of the orthotropic shell, the main axes of 
rigidity of which coincide with the coordinate axes s and 0 , can be 
stated in decoupled form for the section forces Ng, N 0 , Ns 0 for 
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Figure 1. Dimensions and Designations of the 
Truncated Cone 


the section moment\uns Mg, Mq, Ms 0 utilizing the usual assumptions^^'’ 
for the investigation of plates and shells, (In Figure 2 the forces and 
momentiims with positive sign are indicated for an element of the shell.) 



(2. 1) 
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Figure 2. Section Forces and Section Momentum s on an 
Element of a Conical Shell 


Here €g and €0 represent the expansions of the central surface of the shell 
in the direction of the generatrix and in the direction of the circum- 
ference. Ys 0 is the shear of one element of the shell s d0 ds (Figure 
2 ). kg and Vlq are the curvatures; kg 0 is the torsion of the central sur- 
face of the shell. The designates the expansion and shear rigidities, 
the the flexural and torsional rigidities. 

In the further investigations, the first elasticity equations of (2. 1) 
are used in inverse form 




Cii Ci2 0 


"Ns‘ 



Ci2 Czz 0 


Ne 

Vse 


0 0 C 33 

— M 


_Nse 


where the relation between the and the is given by the following 
equations: 
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( 2 . 3 ) 


Cn = 


Cu = 


^22 

Cil Cz2 - Ciz 

-Ciz 

Cu C22 -C12 


C22 = 


C33 - 


^11 

Cn C22 - Ci\ 
1 


'33 


For the isotropic conical shell, the rigidity parameters assume 
the values 


Cu - C22 = ^ , C12 - - ~ 


C 33 =^ 


1 


( 2 . 4 ) 


Du - Dzz - 


Et^ 


12 


(‘-1 


, Di2 - 


vEt^ 


12 (1 - v^) 


2D33 - 


Et' 


12 (l +v) 


where t designates the wall thickness, E the modulus of elasticity and 
G the modulus of rigidity, and v is Poisson's ratio. 

2 . 2 The Kinematic Relations 


The expansions of the central surface of the conical shell 
utilizing the essential members of second order in the displacement 
derivatives® and taking into account the predeformation can be 
expressed by the following relations: 


e Q - 


8u 1 /Bwr 8 wa 

8w 

0s 2 y0sy 0s 

“8^ ’ 

u - w cot a j 1 


s s sin a 

8§ ^ 


1 8w 8 wa 

8 0 


is sin a/ (j y 86 

8v _v ^ 1 8u 

~ 8s s s sin a 80 


+ 


1 8wV 

2 Is sin a d dj 


8w / 1 8w\ 

8s Is sin a 8 


+ 


1 / 8w A Bw ^ 8wa 8w\ 

s sin a I 80 8s 8s 8 6y 


> ( 2 . 5 ) 


Here u, v, w designate the local displacements of a point along 
the generatrix in direction of the circumference and perpendicularly 
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with respect to the central surface of the shell (Figure 1), ( A - 

initial deformation) represents a predeformation of the shell perpen- 
dicular with respect to the central surface of the shell in the same 
sense of direction as the displacement w. 

The displacement-expansion relations of the ideal conical shell 
(wA = 0 ) reduce to the familiar equations of the theory of thin plate s^^ 
for the case a= 90®(si— ►rj, S 2 ~^^Zf ds“^dr). 

Making the transition to the cylinder (a-^0, s— ►oo, s sin a-^R, 
ds— ►dx) from (2. 5), the relations follow which Marguerre^ and Donnell^ 
have derived for the predeformed cylinder shell. 

With Donnell’s simplifications^ extended for the conical shell^, for 
the curvatures of the central surface of the shell follows 

0^W 

ks = 

k - -i. ^ 3^w 

^ 2 “2 
s 0 s s sin a 30 

0 s \s sin a 0 0 / 

By means of the equations of elasticity (2. 1) and the relations (2. 5) 
and ( 2 . 6 ), the section forces Ng, N 0 , Ng 0 and the section momentums 
Mg, M 0 , Mg 0 can be expressed with the displacements u, v, w, and 
their derivatives. The differential equations from which the displace- 
ments u, V, w can be determined can be derived from equilibrium 
considerations on the deformed shell element or by means of the 
calculus of variations from the energy. In the present discussion the 
second way is used. 

2. 3 The Potential Energy of the Deformed Conical Shell 


r (2-6) 


The total potential energy of a conical shell deformed by edge 
loads and inner pressure 


n = III + Hz + 113 + 1I4 + 115 


(2.7) 
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is composed of the following contributions: The expansion energy of 
the central surface of the shell 


02 S2 


~ \ |Ns€s + N 0 € 0 +Nse'yse| s sin a ds d 0 , 


01 Si 


the bending energy 
02 S2 


II2 - — J' J' I Mg kg + M 0 k0 + 2 Mg 0 kg0| s sin a ds d 0 , 

01 Si 

the potential of the edge load at the edges s = constant 
02 


II 3 = - y |s |NsU + Ns 0 
01 ^ ^ 


, 0w - 1 9w 

V + Ms + Ms 0 

3 s s sin a 3 0 


+ 


Qswjj 


S|2 

Si 


sin a d 0 , 


the potential of the edge load at the edges 0 = constant 


H4 



N0 V + Ng0U + M0 


1 3w 

s sin a 3 0 


M 


3w 

s 0 

3 s 


+ Q0W 



Bi 


and of the potential energy of the inner pressure 


02 S2 


n< 


=// 


pw s sin a ds d 0 


6] 


Si 


9 


where Ng . . , Mg. . , Qs • • , Qq are the edge loads. 

2.4 The Differential Equations of the Stress and Deformation State 


2.4. 1 Derivation of the Differential Equations 


The system of differential equations describing the 
state of the stresses and deformation of the orthotropic conical shell 
is derived via the potential energy, the first variation 5 II of which must 
go to zero for all states of equilibrium of the deformed shell. With 
(2. 5), (2. 6), and (2. 7) and the substitution (2. 8), 

0 = 6 sin a , (2. 8) 

the first variation of the total energy is formed and is arranged accord- 
ing to the displacement differentials 6u, 6v, 6w. 

In order to satisfy 6 II = 0, the double integrals over the surface of 
the shell have to vanish individually since 6u, 6v, and 6w are independent 
virtual displacements. This leads to the following equilibrium con- 
ditions: .. 



The boundary conditions at the edges s = constant and 0 = constant 
are obtained from the single integrals in 6IL For the boundary con- 
ditions at the upper ( s = Sj) and lower (s = S 2 ) edge of the shell 
follows: 
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■KH i TVT /SW 9wa\ , T-T /9W . 8wa\ 

Me + sN3(^ + -ar) MTi- + “ar) 


■ J7 (®^s) " ^ ^ ° or 6w = 0 


as 

Ms = Ms 




>( 2 . 10 ) 


Ns = Ns 


Ns0 = Ns0 


or 6u = 0 
or 6v =0 


For the closed conical shell (0^0^ Zir sin a), which is treated here, 
the boundary conditions at the edges 0 = constant are replaced with the 
periodicity requirement with the period 2tt sin a. 

In (2. 10) the general static and geometrical boundary conditions of 
a conical shell are given for arbitrary edge loads. The further 
investigations are restricted to shells under axially symmetric load 
(axial load and inner pressure, Figure 3). 



Figure 3. The Load of the Conical Shell 
Investigated 


1 

Corresponding to the test conditions , the following necessary 
geometrical conditions shall apply to the edges of the shell caused by 
the stiff end-disks bordering the shells at the edges s = si and s - Sz 
(Figure 4): 
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Figure 4. Displacement of the Lower Edge of the Cone 
(s = S 2 ) for Constant Radius of the End 


For such a conical shell under axially symmetric load, the follow- 
ing axially symmetric boundary conditions are obtained at the edges 
s = Si and s = S 2 . 


w = const 


u 

cot a 


Ms = 0 


or 


0W 

a s 


= 0 


N 


s 


= N 


s 


or 


u □ const 

V = 0 


( 1 ) 

(2a) 


(2b) 

(3a) 

(3b) 

(4) 


> ( 2 . 12 ) 


In order to be able to investigate the effect of the boundary con- 
ditions on the buckling load, the condition (2a) --pin- jointed edge--as 
well as the condition (2b)--fully restrained edge--is taken into account 
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in the further considerations. The boundary condition (2b) corresponds 

13 • ^ 

to the test conditions since the shell has been lap- glued to the end 

disks. The system of differential equations (2. 9) with the boundary 

condition (2, 12) generally describes the stress and deformation states 

of the axially loaded orthotropic conical shell. These are: 

1) The state of the shell before buckling- -here called basic state 
of the stresses. 

2) The states of equilibrium at the branching point starting from 
which the shell buckles to assume another configuration- -the edge load 
Ns of the branching point is called critical load or buckling load in the 
further investigations. 

3) The postbuckling range. 

The following considerations are restricted to the determination 
of the basic state of the stresses and to the determination of the branch- 
ing point. For these states (1 and 2) special forms of the system of 
differential equations (2. 9) can be given. 

First the state of the shell prior to buckling shall be investigated. 

2.4. 2 The Differential Equations for the Stress and 
Deformation State Prior to Buckling 


In the further investigations a predeformation w^(s) 
is taken into account which is independent of 0 and includes the special 
case w^ = 0 (ideal cone). 

Considering first the state of the shell for small load (axial load and 
inner pressure) it can be assimied that the stresses and deformations 
are small of first order. Thus the system of equations (2. 9) is linear- 
ized. For the system of linear differential equations with the boundary 
conditions (2. 12) only one solution exists. 'Under the given boundary 
conditions (2. 12) the solution is independent of the circumferential 
angle 6, it is axially symmetric. For the system of differential equa- 
tions of the isotropic conical shell obtained from (2. 9) by linearization, 
a solution in closed form is given in Section 3. 1 of this paper. 

In the following part of the paper, the axially symmetric basic 
state of the stresses shall be investigated for large loads as well (e. g. , 
in the region of the buckling load). Thus the nonlinear members in the 
basic equations- -corresponding to equation (2. 9)--are taken into account. 
Because of the axial symmetry, all derivatives of the section forces and 
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section momentiims and deformations with respect to 6 go to zero. 
Thus the system of differential equations (2. 9) simplifies 



with the boundary conditions (2. 12). The index ( )o of the section forces, 
section momentums, and displacements designates the basic state of 
the stresses. 

Ns0o ^ ® Vo = P (2. 14) 

respectively. 

Due to the axial symmetry of the stresses and' deformations, vq 
and Ng 0 ^ become identically zero. Thus the differential equation 
(2. 13. 2) and the boundary condition (2. 12.4) are satisfied. The two 
differential equations (2. 13. 1) and (2. 13. 3) represent the basic equa- 
tions of the axially symmetric state of the stresses which in the further 
investigations shall be described by the section forces Ng^ and the 
displacement Wq . 

From equation (2. 13. 3), the section force can be eliminated 

by means of (2. 13. 1). After integrating once from (2. 13. 3) follows 

- sNs„(cot +T K„ = 0 , 

(2. 15) 

where Kq designates an integration constant still to be determined. As 
additional differential equation, the compatibility condition from the law 
of elasticity ( 2 . 2 ) is derived by replacing the expansions €so ^ind 
with the displacements uq and wq corresponding to (2. 5) and by then 
eliminating uq . Taking into account the relationship (2. 13. 1) it follows 
that 
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(2. 16) 


-Cz2 




ds^ 


3C22 S 


dNso 

ds 


+ 


^Cii 



Nso 


dwp 

ds 


Icot a 



dwp 

ds 


+ 


dwA\ 

ds f ' 


With (2. 1) and (2. 6\ MgQ and Mqq in (2. 15) are replaced with the 
displacement wp and its derivatives. Then two coupled nonlinear 
differential equations for Ngg and wp are obtained in which wp appears 

/d^w \ 

only in its derivatives with respect to s I ^ ; 9 = 1» 2, 31 and the 

axially symmetric prebuckling wa only in the form dwA/ds. 


To reduce the order of the differential equation (2. 15) it is thus 
substituted: 


dwp 

ds 


= To ; 


dwA _ 


ds 


= qA< 


(2. 17) 


The resulting system of differential equations is further simplified 
by the coordinate transformation 

z 

s = s,e O^z^zp 


zo 



(2. 18) 


A system of two nonlinear coupled differential equations of second order 
for qp and Ng^ is obtained, 

-C 22 N" 3 ^ - 3C22N'g^ + (Ca - C 22 ) Nsp = qp (cot a + y To + 


Du q'i - D 22 qo = sie^ Sie^ NgQ (cota + qp + q_^ -^psfe^^ _ Kq 


(2. 19) 


where ( )' = — - . 

dz 


Introducing the dimensionless functions Np and Qp , the new 
parameters (C 22 , Ca, C 12 ; Dii, D 22 » D 22 ; sj zp, cot a) essential for 
the behavior of the solution of (2. 19) with the boundary conditions (2. 12) 
can be combined to six characteristic quantities {Cn/Czz, C 12 /C 22 , 
D 22 /D 12 , D 12 /D 11 , Zp , X|. With the substitution 
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1 


sie Si 


„ _ C 22 p, 

qo Wo 


„ _ -VDii C22 ^ 

Qa ’ 


system (2. 19) becomes 


N';-|^N„= -e^ 


Qo (x + -^ Qo + 


0 D 


( 2 . 21 ) 


^ Qo = |no (X + Qo + Qa) " + Ko)| , 


Si cot a 


Dn C 


n '->22 


1 

P =— P 


2 cot a 


( 2 . 22 ) 


cot^ a 


The solutions Nq and Qq of the system of differential equations 
(2. 21) must satisfy the boundary conditions (2. 12. 1) and (2. 12. 2). Due 
to the axial symmetry--vo = 0--condition (2. 12. 1)-- wo = Uo/cot a-- 
means that the circiimferential expansion at the edges z = 0 and 
z = zq becomes zero. With the aid of the substitution (2. 17) and (2. 20) 
and the coordinate transformation (2, 18), the boundary conditions 
(2, 12. 1) and (2. 12. 2) are expressed by the variables Nq and Qo : 
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Cz2 N*o + Cj2 No - 0 


( 1 ) 


(this means €0^ = 0) 
1^11 Qq + Di2 Qo = 0 

(pint jointed edge) 


or 


Qo = 0 


(2a) 


> (2.23) 


(2b) 


(fully restrained edge) 

From the conditmn (2. 12. 3a) and (2. 12. 3b) respectively, the 
integration constant Kq in (2.21) is determined. For the range of the 
axially symmetric state of the stresses up to buckling which will be 
investigated in this portion of the paper for each edge load Nq there 
exists a one-to-one correlation of an edge displacement uq . For the 
numerical evaluation Kq is given, and from the solution of the system 
of differential equations (2.21) the edge displacement Uq and the edge 
load No are determined. 

The determining equation for the displacement wq perpendicularly 
with respect to the central surface of the shell is derived from the 
relation (2. 17), 

z 

Wo =y* ybn C 22 Qo (C) dg + Kj . (2.24) 

0 


Without affecting the generality, the integration constant Kx can be set 
equal to zero. This means that the edge z = 0 is retained while the 
shell is deformed, while the edge z = zq can move (Figure 4): 


o 

II 

o 

o 

II 


and 

Uo = 0 

z = zo : Wo = Wo = 

Uq 

and 

Uq = Uo 

cot a 


2.4. 3 The System of Differential Equations for the Additional 
Stresses and Deformations at the Branching Point 

The axially symmetric basic state of the stresses which 
has been described in the last section is maintained as long as the edge 
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load No is smaller than the critical load Ncj.. When reaching the 
critical load (branching point) at least two positions of the equilibrium 
exist, one axially symmetric (identical with the basic state of the 
stresses) and one more or less nonaixially symmetric (starting of 
buckling). At the branching point the smallest external perturbations 
render the basic state of the stresses unstable; a nonaxially symmetric 
state of equilibrium is assumed. 


In this portion of the paper the system of differential equations 
shall be derived from which the nonaxially symmetric stresses and 
deformations occuring for buckling can be determined. 


With the equation of elasticity ( 2 . 1 ), the relations concerning the 
state of deformation ( 2 . 5 ) and ( 2 . 6), and the equilibriiim relations 
( 2 . 9 ), the section forces Ng, N0, Ns0, the section momentums Ms, 
M0, Mg0, the deformations €s, « 0,^80' ^s» ^^0» ^s0, and the dis- 
placements u, V, w can be determined for every state of equilibrium 
of the shell. By introducing a stress function, 


N 


s 


i 11 + J_ 

s 0s s^ 00^ 


N0 


l!l 

ds^ 


Nse =- 


a 


0 s 




^ ( 2 . 26 ) 


which, corresponding to equation (2. 26 ), is defined in such a way that 
the two first equations of (2. 9) are identically satisfied; the system of 
nonlinear partial differential equations can be reduced to two equations 
for the stress function f and the displacement w. The first relation for 
f and w follows from equation ( 2 . 9 . 3 ) by expressing the section momen- 
tums by means of the law of elasticity (2. 1) and by means of the equa- 
tions (2. 6) by w and its derivatives, and the section forces are elimi- 
nated by means of (2. 6). 


Due to the condition of compatibility, the second differential equa- 
tion for f and w is given. From the relations concerning the state of 
deformation ( 2 . 5 ), u and v are eliminated. The expansions €g and €0 
and the shear ^39 are replaced by means of the law of elasticity (2. 2) 
and the substitution ( 2 . 26 ) with the stress function f. Equation ( 2 . 9 . 3 ) 
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and the condition of compatibility form a system of two nonlinear, 
partial differential equations of fourth order in s and 0 for the stress 
function f and the displacement w: 


Lo(w) 


= /i ii + I ili\ , a^wA 


(i il 

\s 0s 


00^ 
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(2. 27) 


L<o(y) represents a linear differential operator: 


0g^ 


i l!l 

8 0s' 


S^ 0s' 




0 y 


03^00^ 


1 Jh\ 

0s0eV 


s^ 0s0e^ 
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where 


19 



af - Czz 


ci^ - Dn 

|3w = D 22 Pf = Cii 

"Yw “ 2 (Di2 + 2 D 33 ) 'Yf = 2Ci2 + C33 . 

The boundary conditions for the system (2. 27) are given in (2. 12). 

The system of differential equations (2. 27) generally describing 
the states of stresses and deformation of the orthotropic conical shell 
agrees with the system of differential equations given by Schnell® for 
the ideal isotropic conical shell (w^ = 0)* 

As described in Section 2. 4. 2, the conical shell first deforms 
axially symmetrically for the given geometrical boundary conditions 
and under the load of axial load and inner pressure. When reaching 
the critical load nonaxially, symmetric contributions are added to the 
axially symmetric stresses and deformations so that for the branching 
point it can be set: 


w = Wo (s) + wi (s , 0) 
f = fo (s) + fi (s, 0) , 


(2. 28) 


where fo and Wo designate the stresses and deformations of the axially 
symmetric basic state of the stresses. fi and wx are the nonaxially 
symmetric stresses and deformations which are added when making the 
transition to the buckled state. If these expressions (2.28) are inserted 
in the system of differential equations (2. 27) it is obtained: 
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and 


9^wi ^ 9^Wo ^ 
9s^ 


I \ / 1 9 fi 1 9fo 1 9^fi\ ( 9^wi 

Lo W 0 + wi) = + + - — T- 

\s 9s s 9s 90 / \9s‘' 
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\2 0s" 
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3s^ / \s^ 30 s 3s00y 


/9^wi ^ 3 ^wq\ 

\3s^ 3sV 


1 / 1 3wi 1 3wo 

— (cot a + + 


2 0s 2 0s 


(2.29) 


9w^\ ^ 1 1 3^wi 




0 s / 2 002 


The axially symmetric contributions f q , wq may be expressed with 
the relations (2. 26) and (2. 17) by NgQ and qo . 

For the branching point, system (2. 19) can be split into two deter- 
mining equations for Ngo and qo (and fo and Wq , respectively) --identical 
with equations (2. 19) of the axially symmetric state of stresses and 
deformations- -and into two determining equations for the nonaxially 
symmetric contributions fi and wi which may be written in the following 
fo rm: 
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(2. 30) 
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If we restrict ourselves to the determination of the buckling loads, 
it can be assumed that in the region of the branching point the functions 
fi (s, 0) and wi (s, 0) as well as their derivatives are small of first 
order. Neglecting members of second order, the following system of 
linear differential equations is obtained for the nonaxially symmetric 
contributions fi and wi. 


Lo(wi) = Bo (wi , s Nso) + Bo (fi , cot a + qo + qA) 

Bo(fi) = - Bo |wi , cot a + qo + qA) 
where Bo |yi, cpo) represents the following differential operator: 


(2. 31) 
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= s |?!i; -i (Po +/i !li + i 1 — • 

( 0 s^ \s 8 s aev s 8 s 

In order to obtain an overall picture of the essential parameters of 
the system, dimensionless functions Fi and Wi are introduced. 

With equations (2. 20) with the substitution 


Wi - e ^Dii Czz Wi 


fi - e Dll Fi 


(2. 32) 


and with the coordinate transformation ( 2 . 18), the system of differential 
equations (2. 32) becomes 


Li (Wi) = Bi (Wi . No ) + Bi (Fi , \ + Qo + Qa) 
Li (Fi) = -Bi(wi, X +Qo + Qa) . 


(2. 33) 


where 


F-i (yi) = ^ yi +— yi + — ^yi + 2yi + yT^ 

ay ay ay 

is a differential operator with the constant coefficients 
<IF “^22 > ' "YF ~ 2 Ci 2 + C 33 

a-yy = Dll , p-yy = D 22 , "YW" ^ (D 12 + 2 D 33 ) 

and 


/ \ Z 


<u 

II 

0 

e- 

cpo yi + ^<Po + cpo'jyi' + cpo' \y + yi / 


In (2. 33) the differential quotients 8 ( )/0z and 8 ( )/00 have been 
replaced with the following symbols: 


ill 

0Z 


( )' 


and 


ill 

00 


0 . 
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The determination of the nonaxially symmetric stress and defor- 
mation and contributions Fi and Wi is performed under the condition 
that in the instance of buckling- -during the transition from the axially 
symmetric basic state of the stresses to an infinitesimally adjacent 
nonaxially symmetric state of the equilibrium- -the axial edge load Ns 
(2, 12. 3a) and the radial edge displacement w (2. 12. 1) remain constant. 
The axially symmetric stress and deformation contributions No and Qo 
satisfy the boundary conditions corresponding to (2. 12). Thus the 
additionally occurring nonaxially symmetric stresses and deformations 
appearing in the instance of buckling must satisfy the following homo- 
geneous boundary conditions at the edges z = 0 and z = zq: 


Wi = 0 

(1) 

0 

II 

(2a) 

Wi' = 0 

(2b) 

(Ns)i = 0 

(3) 

Vi H 0 

(4) 


> (2.34) 


Index { )i designates the additional stresses and deformations 
occurring in the instance of buckling. 

From the boundary condition (2. 34. 1)--Wi = 0, i. e, , also 3Wi/30 
= 0--and condition (2. 34.4)--vi = 0, likewise 3vi/30 = 0 - -it follows, due 
to the necessary geometrical condition (2. 17), that for the edges z = 0 
and z = Zq 

(^ 0)1 = 0 because (^s)i = 0 follows that ^Ne)i = 0. (2. 35) 

With the aid of the law of elasticity (2. 1), the definition of the stress 
function (2.26), and the substitution (2. 32) from (2. 34) and (2. 35), the 
following boundary conditions for functions F] and Wi can be formulated: 
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11 


Dll + Di 2 


or 


Wi 

= 0 

(1) 

Wi" + Wi' + wi' 

= 0 

v2a) 

Dll + Di2 


means 

Wi' 

= 0 

(2b) 

Fi' + Fi + Fi" 

= 0 

(3) 



means 

Fi" + Fi' 

= 0 

(4) 


(Ms)i = 


>(2. 36) 


means 


(Ns)i = 


(Ne). = 0 


Since the system of linear partial differential equations (2. 33) as 
well as the boundary conditions (2. 36) are homogeneous, the deter- 
mination of the nonaxially symmetric stress and deformation contri- 
butions Fj and Wj is an eigen- value problem. Fj and Wi are the 
associated eigen functions; the axial edge load No (and the edge dis- 
placement uo or the constant Kq , respectively) of the axially symmetric 
basic state of the stresses represents the eigen-value of the system. 


To solve the system of differential equations (2. 33), the functions 
Fi and Wj which are periodic with 2n sin a in the direction of the 
circiomference are expanded with r| = n/sina into series of the form 


00 

Fl = ^ Fin (z) cos T)0 
n=i 


00 


Wi =23 
n=i 


Win (z) cos T| 0 


(2.37) 


and are inserted into the system of differential equations (2. 33). By 
comparing the coefficients, systems of linear, homogeneous ordinary 
differential equations are obtained, 
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(2. 38) 


I ->2 + B 2 ^Fijj , X + Qo + Qa^ 

Lz (f m) — ^2 (w,n . X + Qo + Qa) 

where Bz ^yin . cPo^ represent the ordinary differential 

operators, 

X / \ ( 6 ) ay + Py + 6 y " ^ Py A 2)2 

L>2 |yinj= yin' ' yin + / yin 

^ /uTr 


with ap = C 22 . Pf = Cii , 6p = ^2 C 12 + 

o-y^ = Dll , p-^ = D22 1 6w ~ ^ (di2 + 2D33J 
Bz^yin. fPo^= e cpo yin + (<Po + cpo) yin + cPo (l - t|^) yin 


The functions Fm (z) and Wm (z) have to satisfy the following 
boundary conditions according to (2. 36): 



0 

II 

CJ 

(1) 


Dll 

w'ln + Win = 0 

(2a) 


Dll + Di2 


or 

5 

II 

0 

(2b) 

^ (2.39) 

F'in + 

(1 - Fin = 0 

(3) 



F in + F in = 0 • 

(4) 


For each n (number of bucklings in the direction of the circum- 
ference in the instance of buckling) a system of differential equations 
(2. 38) with the boundary conditions (2. 39) exists, the solutions of 
which Fin and Wm represent the eigen functions with the associated 
eigen-value No (n). 
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3. Solutions and Solution Methods 


3. 1 Exact Solutions of the Simplified System of Differential 

Equations for the Axially Symmetric State of Stresses and 
Deformation 


The determination of the buckling stress, i. e. , a solution of 
the system of differential equations (2.38), assvimes the knowledge of 
the axially symmetric stresses and deformations (No, Qo) of the 
unbuckled region--here called the basic state of stresses--corresponding 
to equation (2. 21). For the cylinder shell^, the differential equations 
of the axially symmetric state of stresses and deformations can be 
solved linearly and in closed form. The corresponding system of 
differential equations of the conical shell ( 2 . 21 ) remains nonlinear 
even in the axially symmetric case. A general closed solution of this 
nonlinear system of differential equations for the basic state of the 
stresses is not possible with the known methods. Therefore, first, 
two methods shall be stated (a and b) by means of which solutions for 
No and Qo can be found when the system of differential equations is 
simplified. 

a) Diaphragm stress state 

Under the assumption that the deformations of the shell 
perpendicular with respect to the central surface of the shell are con- 
stant over the entire shell of the cone--this means that Qo is identically 
zero- -and that no prebuckling occurs (Qa = 0)», No is determined from 
equation ( 2 . 21 . 2 ): 

No = x(pe^^ + Ko) . (3.1) 

With the abbreviations (2. 20) and (2. 22), this becomes 



2 cot a s cot a 


the solution for the flexural-free state of diaphragm stress correspond- 
ing to Schnell's^ Equation (2) and Schiffner's® Equation (7) if the axial 
load (Figure 3) is set 

(3.3) 


p = -2 tr sin^ a Ko . 

For N 0 p follows then from equation (2. 13. 1): 
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(3.4) 


N 


00 “ 


p.s 
cot a 


The diaphragm stresses Ns^ and N 0 ^ corresponding to equations 
(3. 2) and (3.4) satisfy the equilibrium conditions (2, 13) and the 
boundary conditions 

Ng^ = Ng at the edges z = 0 and z = zq , 

not, however, the compatibility condition and the boundary condition 
^00 “ 0 thus are no valid solution of the boundary value problem 
investigated here. 

b) Linearization of the System of Differential Equations 


For small edge loads and small inner pressure, the 
deformations and stresses of the ideal cone (Qa ~ small of first 

order. If, therefore, the system of differential equations (2.21) is 
linearized- -this linearization means a neglect of Qq compared with \-- 
the following system of differential equations is obtained: 


No“ - ^ No = -e^ Qo\ 


*21 


Qo” - ^Qo = \ 


(P 


No - X(p +,Ko) 


(3.5) 


and after eliminating Qq, respectively 


r I 


No 






No 




0 


= - e 


2 z 


Nf 


\(pe^^ + Ko) 


(3.6) 


Equations (3.5) and (3.6), respectively, correspond to the dif- 
ferential equations of the linear bending theory^^ for axially symmetric 
stresses and deformations (container theory). 


For the isotropic conical shell the rigidity parameters of which 
may assume the value Cu/Ci 2 = D 22 /D 11 = 1 according to (2. 4) the dif- 
ferential equation (3. 6 ) can be solved in closed form. As for the 
cylinder shell, the particular solution of ( 3 . 6 ) is identical with the 
diaphragm stress (3. 1 ). Together with the solution (modified cylinder 
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function- -Kelvin function^^) of the homogeneous differential equation 
it is obtained 


0 - ^ ^ ber 2 ^ 2 \/\s ^ + Kj beiz^Zv^t^^ 

+ K 3 her 2 ^ 2y/^^ + K 4 hei 2 ^ 2 v/\s ) 


(3.7) 


For the range of parameters (Section 4. 1 ) investigated here-- 
argument x = 2 ^s > 10 ^, order “ = 2 --these functions are not tabulated. 
However, J ahnke give s a first approximation of the Kelvin functions 
for the parameters x » 1, x »T7 . Inserting this approximation into 
solution (3. 7) for No is obtained 


No = X ( s^ + Ko W ^ 


sf 




sin h/^ 


Ki sin Jz\s 


+ Kz cos y/zXs 


[f cosh y 2 \s 


K 3 sin \/z\s + K4 cos y/z\^ 


(3.8) 


Here the constants , Kz , K 3 , K 4 are to be determined from the 
boundary conditions (2. 23) and Kq from the condition (2, 12. 3a) or 
(2, 12, 3b). For the numerical evaluation, however, Kq is given so that 
then the edge load Rq sind the edge displacement uo is calculated, 
respectively. 


The approximate solution (3, 7) and (3. 8 ), respectively, of the 
linearized system of differential equations from ( 2 . 21 ) represents a 
good approximation for the basic state of the stresses for small 
edge loads and small inner pressure. However, it cannot be used for 
determining the buckling load since the system of differential equa- 
tions (2. 31) for the additional stresses and deformations fi and Wi 
appearing in the instance of buckling has been derived under the con- 
dition that No and Qo (and Ng^ and wo , respectively) are solutions of 
the nonlinear system of differential equations ( 2 . 21 ). 


3. 2 Approximate Solution of the Linear Buck ling Equations by 
Means of the Energy Method 


In the following portion of the discussion, the state of the 
diaphragm stress (3. 1 ) described in Section 3. 1 shall be applied to 
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the determination of the buckling load. Inserting the diaphragm stresses 
(3. 1 ) into the system of differential equations (2. 38), the buckling equa- 
tions of the linear theory of shells is obtained; 


Le 

(w.„) 

= Bz (Win, No) 

+ B 2 

(f in > x) 



Lz 

(f„) 

= - Bz (Win.x) 




^ (3.9) 




z 

/II 1 \ 



where 


Bz (Ym, X) = 

Xe 

(Yin + Ymj 



and 


No = X (p 

+ Ko' 
> 

). 




The system of differential equations agrees with the buckling equa- 
tions given in several source provided the appropriate sub- 

stitutions are observed. 

Except for the special case r\ ~ 0 (ring buckling) which has been 
treated by Seide^^ for the isotropic conical shell with the aid of further 
simplifications, no solution in closed form is known for the system of 
differential equations (3. 9). 

It is stated once more that the state of the diaphragm stress does 
not represent an axially symmetric solution of the shell problem 
investigated here since the compatibility condition as well as the 
boundary condition € 0 ^ = 0 are not satisfied. The further investigation 
according to the linear theory has only the aim to derive a reference 
parameter for the buckling loads calculated according to the nonlinear 
shell theory. 

Since a general solution of the linear buckling differential equa- 
tions (3.9) is not known, an approximate solution according to the 
energy method is determined in the following. An appropriate formu- 
lation for the deformations perpendicular with respect to the central 
surface of the shell is chosen, and from the condition of compatibility 
a stress function is determined with the formulation for Wxn- The 
free values Ak of the formulation are determined from the condition 
for the energy minimum 0 ir/ 0 Ak = 0 corresponding to the variation 
requirement 6 tt = 0 . 

For a basic state of the stresses according to equations (3. 2) and 
(3. 4) --diaphragm stresses- -the following expression is obtained for 
the total energy of the deformed shell in the branching point. 
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(3. 10) 


where the first two double integrals represent the deformation energy 
of the forces and deformations occurring additionally during buckling. 
Since these additional forces and deformations are small of first order, 
the nonlinear contributions in €g^, ysQi neglected 


(« 


(L) 

Si » 



(L) 


)■ 


For the displacement Wm appearing in the instance of buckling 

perpendicular with respect to the central surface of the shell, the 
formulation 


Wjn = Ae ^ sin k z with k = ^ - - (3. 11) 

zo 

is chosen. This formulation satisfies the boundary condition (2. 39. 1) 
at the edges z = 0 and z = zq ; Wm = 0 and corresponds to the so-called 
’’checkerboard buckling configuration” with ffi-half-wave,s in the longi- 
tudinal direction and n-waves in the circiimferential direction, where 
because of s = sj e^ the half-wavelength increases along the coordinate 
s. For a given the stress function can be found from the compati- 
bility condition(3,9.2) by integration. It is composed of aparticular solution, 


31 



and of the solutions 



of the homogeneous differential equation. 

With these solutions the boundary conditions for the stress function 
can be satisfied. Investigations on isotropic conical shells, however, 
indicated that the effect of the boundary conditions on the buckling loads 
is small according to the linear theory. For this reason only the 
particular solution is taken into account. 

Inserting (3. 11) and (3. 12) into (3. 10) and after integration, the 
total energy of the orthotropic conical shell is calculated from (3. 10), 

From the minimum condition 

(3. 13) 

9 

in addition to the trivial solution k = 0 (no buckling), 
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1 + 2 k‘ 


K'-n 


+ 4 Tl' 


is obtained for the force at the edge z = 0 (tension positive according 
to Figure 2). 

A general evaluation of (3. 14) is involved while by restricting to 
certain shell parameters equation (3. 14) can be essentially simplified. 
In the following only shells shall be considered the geometry of which 
satisfies the equation 


zo = i n ^ ^1 

Si 


S 2 

Si 


^ e 


(3.15) 


and the opening angle a of which is smaller than 45® (see Figure 5). 
For this range of the shell parameters it is 
k^ » 1 , rf» 1 
so that it may be set 
k^ ± 1 « k^ 

± 1 * ti" . 
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Figure 5. Possible Cone Configurations With 
1*2 /ri = S 2 /si = 2 and e, Respectively, for 
Various Opening Angles of the Cone a<45 ° 


From (3. 14) follows then 
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So 


-Zo 


Zo .2 
— cot a 


T|^ C n + C 11 + C 22 + ^2Ci2 + C33^k t| 



If first the special case of the axially symmetric buckling con- 
figuration (ring buckling) is investigated, i. e. , t| = 0, from 


0N 

0k 


Sfi = 0 
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. 16 ) 
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a number of the half-waves in longitudinal direction can be determined 


= / ■ ^ 6 j cot a 

VC22 Du (1 -e"'‘^°) 


Cii 

C 22 


with the associated ring buckling stress 


(3. 17) 



+ . L . l) 111. . (3.18) 

2 si ' ^ \ CzzDzz/ 2 cot a 


From this formula it follows that the buckling load increases 
linearly with the inner pressure. This solution agrees with the ring 
buckling stress of the linear theory of stability stated by Schnell^ for 
the special case of the isotropic shell. 

The general investigation of the buckling equation (3, I 6 ) is 
restricted to isotropic and longitudinally reinforced orthotropic shells. 
For these shells it is 

Cii^Czz; Dii>D22 (3.19) 

and because of (3, 15) 

Cii + C 22 « C 22 k^ 


Djj k^ + D22 ~ k^ . 


For a further calculation, it turns out to be practical to combine 
the orthotropic rigidities to three rigidity characteristics'^’^^’^^ : 
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Ciz + 2 C 33 


\/^Trc 


i?p= 
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Di2 + 2D23 
\Dii D22 
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11 '-'IX 


(3. 20) 


D 22 C 22 


= disc characteristic, = plate characteristic = main rigidity- 
characteristic. For the isotropic shell it is ~ *^p ~ "YN ~ Equa- 
tion (3, 16) can be further simplified with (3, 19) and (3. 20): 
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where 
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c|)i = 1 + 2 ,Js 


Xo 




(3.21) 


Now the minimiam buckling load is determined in dependence of the 
wave number ratios Xo and |ji. From 


3N 


So 


axo 


= 0 


follows 


Xo = ^ J cot 


V 2zo 

TTT^ 


C 22 D 22 4>2 


(3.22) 


In (3. 23) only the positive root is of significance since Xo < 0 
results in complex m and n which is geometrically meaningless. 
With (3. 22) it follows then that for the critical buckling load 
(dimensionless load parameter) 
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cr ~ 



+ P 



(3. 23) 


where 
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For the conical shells (zq ^ 1) considered here, the value of the 
root dr is about 1, 


Zo 

1 


2 

1 

-Zo 

- e 


« 1 , 


and when making the transition to the cylinder (zq 0) it is exactly 1, 

The ratio of the wave numbers |j. is determined from the minimirm 
condition 0(rci./3(jL^ = 0: 


(3.24) 
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For a given orthotropic conical shell satisfying the geometrical 
requirements according to (3. 15) and the rigidities of wbich are within 
the limits given by (3. 19), from (3. 23) the critical buckling stress o-^r 
can be determined for any inner pressure where |jl^ is to be determined 
from equation (3. 24). 

Observing the parameters matched to the geometry of the cone 
equations, (3, 23) and (3. 24) are identical with the buckling conditions 
of the longitudinally compressed orthotropic cylinder shell discussed 
by Schnell^^ Therefore equations (3. 23) and (3. 24) are not further 
investigated in this paper. 
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3, 3 Iterative Solution by Means of the Difference Method 


Closed solutions of the system of nonlinear differential equa- 
tions (2.21) for the basic state of stresses and of the eigen-value problem 
(2. 38) for the stress and deformation contributions in the instance of 
buckling are not known. Thus an approximation method for solving both 
problems (2. 21) and (2, 38) is described in this section of the paper. 

The derivatives with respect to z and the differential operators L{ym) 
and B(yin. cpo) --generally described by 

^ /h^ ( 4 ) h^ " h ' \ 

G yn . Vn > Yn . Ynl 

\4! 2! 1! / 

--appearing in the systems of differential equations (2. 21) and (2. 38) 
are expanded into difference expressions of the form 

^ Uk)+2 akUk + 0(h«) = 0. (3.25) 

\ • • • / 

The Uk are the function values at m equidistant intermediate points 
in the interval 0 <zo. The distance h between the points is 


h = . 

m + 1 

The residual member 0(h®) gives the order of the first nonvanishing 
member of G if the Taylor expansions are balanced, i. e. , balancing has 
been performed up to the members with h^ inclusively. We first con- 
sider the axially symmetric state of the stresses and deformation. 

3. 3. 1 Approximate Solution for the Basic State of Stresses 
Finite expressions of the form 
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G U- + 6u Ui)+2 aik Uk + r,o (jj k^Uo +Uo) = 0 

k=l 


for i = 1,2,3 


/2m \ 

G (fr U- + 6u Ui)+2 aik Uk = 0 

k=i-3 


for i = 4, S. . . (m-3) 


>(3.26) 


/h^ " \ > /h ' 

G (jr Ui + 5u Uij aik Uk + ri„,+, k^U^^, 


k=m-6 


+ U. 


m+ij - 0 


for i = m-2, m-1, m 


are derived for the system of differential equations, Ui are the function 
values No and Qo at the m intermediate points; 6u follows from the 
coefficients of the differential equation (2.21), 



6 


Q ~ 



Dz2 

Du 


The boundary conditions (2. 23) are taken into account in the finite 
expressions for the edge-near points (i = 1 , 2, 3; i =m -2, m - 1 , m) in 
the form 


'u — ® 


The factor means 


(i = 0, i = m + 1). 


k 


1 c 

N= h- C 


11 

12 


corresponding to (2. 23. 1) 




h D 


12 


- corresponding to (2. 23. 2a) 
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or 


kQ = 0 


- corresponding to (2, 23, 2b). 


If the differential equations at any intermediate point are replaced 

with the finite expressions formed according to equations (3. 26) a 

system of differential equations is obtained which can be combined to 
• 1 7 

matrix equations : 

(afi) (Qoi) - (bkk) (Noi) - (Qoi Noi) + (ci) = 0 
(alk) (Noi) + (bkk) (Ooi) + ^ ° 


(3.28) 


1 « i, k m 


^a^) represent quadratic matrices 




ail 

_(j) 

a-12 

a-13 • • 

3-21 

Jj) 

a-22 

az3 • • 

a3i 

a .32 

a-33 • • 


with the elements a^j^ which are composed of the coefficients of the 
finite expressions (3.25). 

(Qoi), (Noi), (QoiNoi). (Ooi). (c are vectors the elements Qoi 
and Noi which are the function values Qo and No at the m intermediate 
points and the elements ci of which are determined from the equation 

Ci = (p e + Ko) . 

(bkk) is a diagonal matrix 
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0 



t>22 


t>33 



with the elements 

^kk = ^ + ^Ak • 

Furthermore, the vector can be replaced with 



Qoj 

Noi' 


0 


No/ 

^Qoi Noi^ - 

Q02 

No 3 

“ (Ookkj (^oi) ” 

Q022 


N02 


Q03 

N03 


Q033 


N03 



• 


0 


• 


whereby the following relation exists between the elements Qoj of the 
vector and the elements Qojj the diagonal matrix 

Qoj = Qojj • 


From equation (3. 28. 1) vector (Noi) 'can be eliminated with (3. 28. 2) 
so that a matrix equation^ (Qoi) is obtained for the vector (Qoi). 


This nonlinear matrix equation is solved by means of Newton's 
iteration: For a vth approximation (Qoi^^^)> the matrix equation^ 

is satisfied except for the defects 



(3.29) 


By means of corrections (a Qoi^^ / which are determined from the 
system of linear equations 
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(3. 30) 


8 i (Qo/*'*) 

0Qoj 


A Q,j + 



= 0 


and are added to the vth approximation ^Qoi ^ j, the nonlinear matrix 
equation (3. 28) is iteratively solved step-by-step. The iteration is 
carried out until two successive approximate solutions differ by more 
than a previously specified amount. In matrix form the following 
iteration scheme can be given: 


(4^ (QoI<’'>) 4 (dkk'i (ei'i ^ (Ci) = (Ai'-') 


|(4i) (Vi (aik't ‘ (Vi + (ekk'li (AQoi'i 

= -(A,<i 


M3. 31) 


(Ooi'ik (AQoi'i 


where 

(dk4i = [(•’kk) + (Oofck '"*)j 

(a/i = (Vr[(V(ooi'iV(«»i'^’“)] • 


(o) 

If the initial approximation of this iteration is set Qoi = 0, the 
approximate solution represents an approximate solution of the 

linearized differential equations for the axially symmetric state of 
stress and deformation according to equation (3. 5). 


3.3.2 Approximate Solution for the Eigen-Value Problem 


The differential equations (2. 38) for the determination 
of the stresses and deformations in the instance of buckling are replaced 
with difference expressions of the form given in (3. 25): 
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for i = 1, 2, 3, 4 


i+4 

G|L(ui). B(Ui aikUk = 0 


k=i-4 


for i = 5, 6 . . . m = 4 


>(3.32) 


G jL(Ui), B(Ui,^oi)j +] 


h 

V j ^ik + 3^m+i |9iu 

k“-5 I 1! 


1 f xr” 

Um+ij + Sijn+i p2u Ujn+i 

+ — = 0 
1! J 

for i = m- 3,m-2, 
ra - 1 , m 

The are the function values Fm and Wm at m equidistant 
intermediate points. For <Poi the stresses and deformations of the 
basic state of stresses have to be inserted. 

In (3. 32) the homogeneous boundary conditions (2. 39) were com- 
bined to the form 


92u — u'k + Uk = 0 
1 ! 

h^ " h ' 

— Uk + Uk = 0. 

2! 1 ! 


(3. 33) 
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For the coefficients 9i^, 


92 u it is then 


9iw = 0 (1) 

meaning Wm = 0 (2.39.1) 

92w=— — (2a) 

h Dll 1^22 

meaning Wm + Win = 0 (2.39.2a) 

Dll + Dz2 


or 


5zF = 0 

I 

meaning Wm = 0 


h 1 - T)^ 

meaning F m + (l 



h 



Fin - 0 


(2b) 


(2.39.2b) 


(3) 


(2.39.3) 


(4) 


I II 

meaning Fm + Dm = 0- 


(2.39.4) 


From the difference equation at each intermediate point, 
of linear equations for the F mi and W mi is obtained which is 
in matrix equations. 


(b*iL) (Wmi) = (c?^) (F,ni) - (c^l) (w,„i) 
(bfj (Fmi) = - (clk) (w.ni) 


1 i , k ^ m 


) (3.34) 


a system 
combined 


(3. 35) 


where 



The elements and of the matrices are composed of the 

coefficients of the finite expressions. Nontrivial solutions of this 
system of linear homogeneous equations (3. 35) exist only if the coef- 
ficient determinant vanishes. 


det 



(3. 36) 


The critical eigen- value is the smallest edge load No with the 
associated displacement of the edge uq of the axially symmetric state 
of stresses and deformation for which the determinant (3. 36) become 
zero. Of the eigen- values determined for various niombers n of 
bucklings in the direction of the circvomference the absolute smallest 
yields the buckling values (buckling load Nocr)* 


With the approximation methods given here, the boundary conditions 
can be completely taken into account without difficulty. 

For the numerical evaluation, the nxrmber m of the equidistant 
intermediate points- -i, e, , the number of the unknowns- -is increased 
as long as the difference of the critical values of two successive 
approximations Nocr (n, m) - Nocr (n, m+1) stays above a given error 
limit. 


4, Numerical Evaluation and Discussion of the Results 


4, 1 Dimensions of the Shell and Parameters 


The systems of differential equations (2,21)--for the deter- 
mination of the basic state of stresses--and (2. 38)--for the determi- 
nation of the critical load- -with the associated boundary conditions 
(2, 23) and (2. 39) contain the following eight parameters 


Si cot a Cii D22 2 C12 + C33 2{Di2 + 2D33) 

■ f , ■ — ■■ ■ > " 

VDii C22 C22 Dit C22 Dll 


and 

C 22 


Di2 

Du 


which are composed of the three geometrical characteristics Sj, zq , 
and cot a and the eight rigidity parameters C 22 . Cxi, C 12, C 33, Du, 
D22, Di 2, andD 23 . 
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For the isotropic conical shell, the rigidity parameters of which 
are given in equations (2.4), the eight parameters reduce to 

Si cot a rm TT 

zo, X - ^12 (1 - v2) , V . 

t 

The mentioned shell parameter X is identical with the cone char- 
acteristic stated by Weingarten et al^^ . As can be seen from the 
systems of differential equations (2. 21) and (2. 38) with the boundary 
conditions (2.23) and (2.39) X is not sufficient, however, for the deter- 
mination of the buckling load of the isotropic cone (as stated by 
Weingarten et al^^) but--aside from Poisson^s ratio v--also the ratio 
of the end radii(r 2 /ri = sz/s\ - e^f has an effect. 

The numerical evaluation is performed for conical shells for the 
stability behavior of which experimental results are available ’ . 
Because of the extent of the computation- -compared with the computer 
IBM 1620/23 available- -the investigation was limited to conical shells 
(Figure 6) with constant height hj^ and constant radius tz of the base 
while the half angle a was varied (a-= 10®, 20®, 30®). The investigated 
conical shells have a wall thickness of t = 0. 255 mm. The elasticity 
modulus has been regarded as constant, E = 525 kg/mm^. Poisson's 
number was v = 0. 3. 



Figure 6. Dimensions of the Conical Shells Investigated 
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In Table I the dimensions of the investigated shells are listed. 

For the isotropic shells the shell parameters zo-and \ in the last two 
columns suffice for the determination of the stress and stability- 
behavior. For comparison, the cylinder (a = 0) with radius tz and 
height hit for which the buckling load has been determined by Fischer* 
is used. 


Table I. List of the Dimensions and Parameters of the Cone 


Desig- 

nation 

Cone Geometry 

Cone Parameter 

Si [mm] 

S2 

Si 

cot a 

Zo 

X 

JOO 

00 

1 

00 

0 

00 

JlO 

9. 75877*10^ 

1. 18022 

5. 67128 

0. 165705 

7. 17212-10^ 

J20 

4. 00440-10^ 

1.46030 

2.74748 

0. 378639 

1.42575-10'* 

J30 

2. 00000- 10^ 

2. 00000 

1. 732051 

0. 693147 

4.48913-10* 


In order to be able to compare the results of the calculation for 
various shells, dimensionless load and inner pressure parameters 
(Nv and P) are used in the graphic presentations. The computed edge 
loads at the edge of the cone z = 0 are related to the critical ring 
buckling load of an isotropic cone under pure axial load which had been 
obtained according to the linear theory according to equation (3. 18). 


Nso (p = 0) = - 


2 cot a /zo + 1 //Du\ 

“77" VT 7^ i 


(4. 1) 


sotrope 


For the orthotropic shells, the ring buckling load of the comparison 
shell (J 30, shell of the same weight and of the same geometrical 
dimensions) is used as reference load. 

For the investigated conical shells (0 S zq ^ J?n2) the root can be 

set 





Zo 


+ 1 
- 1 


1 


(4.2) 
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so that 


Nso(P = 0) = -^^^ 

sx f \C 22 / isotrope 

applies for the ring buckling load (linear theory) of an isotropic conical 
shell^®. 

In Table II the exact values of root (4. 2) are listed for the conical 
shells of Table I. 


cot a 


Et" 


ii ^3(1- v^) 


(4.3) 


Table J1 


De sig- 
nation 

Zo 

Jzo + 1 

iT e^o - 1 

JOO 

0 

1 

JIO 

0. 165705 

1. 00114 

J20 

0. 578639 

1. 00594 

J30 

0. 693147 

1. 01967 


Going in the limit to the cylinder (a 0, cot a/si — ► 1 /R, zq 0) 
from the equation (4, 3) follows the so-called classic buckling load of 
an infinitely long cylinder : 


N 


cl 


Et^ 

R 


1 

V3(l - v2) 


0. 605 



R 


(4.4) 


For the dimensionless load parameter Nv> the reference value 


Nv = 


N. 


N 


Sc?(P = 0) 


2X 


(4.5) 


is obtained with the classic ring buckling load according to formula 
(4, 3) and for the orthotropic shells 
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As inner pressure parameter p is used according to formula 

( 2 . 22 ), 

2 \cot a/ f Dll 

In the case of the orthotropic shell for Du, C 22 the rigidity- 
parameter of the isotropic comparison shell are to be inserted, 

4, 2 Axially Symmetric Basic State of Stresses and Buckling 
Loads of Ideal Isotropic Conical Shells 


The axially symmetric basic state of stresses of the isotropic 
conical shell can be described with a closed solution of the linearized 
system of differential equations (3. 5) for small edge load and small 
inner pressure. From Figure 7 where the edge load Nv of an isotropic 
cone J30 is depicted for pure axial load (p 0) as a function of the 
radial displacement of the edge wr (proportional to the compression 
of the cone) and the axially symmetric configurations of buckling along 
the generatrix for various edge loads, it is seen that for small edge 
loads Nv the curves A (linearized, system of differential equations 
(3. 5)) and B (nonlinear, system of differential equations (2. 21)) are in 
good agreement. With increasing edge load, ring buckling appears in 
the edge zones of the cone --analogously to the cylinder shell--which 
can be comprehended only by the solution of the nonlinear differential 
equations (2. 21). 

In Figure 8 curves Nv vs. |Qomax | are plotted for various inner 
pressures for the cone J30 (pin jointed edge) (for comparison the curve 
calculated from the system of linearized differential equations for 
p = 0 is included). First |Qomax| increases linearly with the axial 
edge load, for higher edge load it increases nonlinearly. When reach- 
ing the critical load, the axially symmetric basic state of stress 
becomes unstable (for this reason the curve above the critical load is 
dashed). The shell starts to assume a new state of equilibrium. In 
the framework of the calculation performed here (linearization at the 
branching point) the branching point with the number, n, of bucklings 
over the circumference appearing in the beginning of buckling is 
determined. Concerning the further course of the curves Nv (i« omaxj) > 
nothing can be stated. In the experiment as well as in the calculation 
the decrease of the number of bucklings, n, with increasing inner 
pressure could be observed (however, in the experiment the number of 
bucklings was determined only for the stable postbuckling region). 
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(A) Basic state of stress calculated from the linearized 
system of differential equations (3*5) 
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Figure 7, Axially Symmetric Basic State of Stress of an Isotropic Cone 
J30 Under Pure Axial Load and Pin- Jointed Edge 




— LINEAR, SYSTEM OF DIFFERENTIAL EQUATIONS(3.5) 



of an Isotropic Pin-Jointed Conical Shell for Various Inner Pressures 
(n - Nvimber of Buckling s in Circvimferential Direction) 





The axially symmetric configurations of the bucklings for various 
inner pressures are plotted in Figure 9 for a constant edge load which 
equals the buckling load of the conical shell for p = 0, i. e. , Nv = Nvcr 
at (p = 0). For a pure axial load (p = 0), axially symmetric buckling 
configurations. begin to form beginning at the edges, decaying towards 
the center of the cone. Due to the inner pressure, the radial stress 
increases with increasing radius, i. e. , with increasing z for which 
reason the radial deformation in the edge zone z = zo is suppressed. 

The tendency that with increasing inner pressure the axially symmetric 
buckling configurations are formed only at the edge z = 0 and there 
with increasing amplitude, typically for cones, is seen. For high 
inner pressure, ring buckling could be observed at the edge z = 0 of 
the test shells before reaching the critical load. 

A comparison of the measured buckling loads (four test shells J30 
of the same dimensions'^) shown in Figure 10 with the critical loads 
calculated according to the nonlinear theory for a fully restrained ideal 
cone (curve III) shows good agreement for large inner pressure. How- 
ever, for small inner pressure- -especially for p = 0--the experimental 
values deviate from the theoretical buckling loads for an ideal conical 
shell. The discrepancy between theory and experiment can be essen- 
tially traced to the existence of prebuckling in the test shells. An 
investigation of the effect of the predeformation on the buckling loads 
is performed in Section 4.4 of this paper. 

The theoretical buckling loads of the pin-jointed cone (curve II) 
are- -like for other stability problems (buckling of a rod, cylinder 
buckling) - -below those calculated buckling values for the fully restrained 
shell (curve m), Curve (I) contains the critical loads which were cal- 
culated according to the linear theory for one-member approximation 
formulation (ring buckling) --corresponding to formula (3. 18). Com- 
pared with curves (II) and (III), the buckling loads determined accord- 
ing to the linear theory represent only a course approximation to the 
stability problem investigated here. 

In Figure 11 the buckling loads according to the linear (curve I) 
and according to the nonlinear theory (curve II) for pin-jointed conical 
shells of the saime height and the Scime base radius rz --according to 
Table I--are plotted versus the opening angle a. It is seen, as shown 
by the corresponding experiments^^, that the thus normalized buckling 
load Nyr--equation (4. 4) --is practically independent of the cone angle 
which on the other side means that the linear buckling load (3. 18) used 
as reference quantity properly represents the effect of the cone angle 
on the buckling load. The values for the cylinder shell (curve II, 
a = 0) entered into the diagreim are taken from Fischer^ and 
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Figure 9. Axially Symmetric Buckling Configurations of an Isotropic Conical 
Shell (J30) = Nycr (p ~ 0) for Various Inner Pressures 







(I) Buckling loads according to the linear theory 
[formula (3. 18)^] . 

(II) Buckling loads according to the nonlinear 
theory [formulas (2.21) and (2. 38)] for pin- 
jointed edge buckling load of cylinder JOO^. 

n: Number of bucklings in circumferential 

dire.ction. 

Figure 11. The Effect of the Opening Angle of the Cone 
on the Buckling Load of Pin- Jointed Isotropic Con- 
ical Shells of the Same Height and the Same Base 
Radius 


satisfactorily supplement the here performed cone investigations 
(a = 10®, 20®, 30®) concerning the buckling load as well as with 
respect to the momber of bucklings n. 
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4. 3 Investigation of an Orthotropic Conical Shell 


With the approximation methods given here (Section 3.4) 
orthotropic conical shells can be investigated as well without additional 
difficulties, provided the neutral plane coincides with the central sur- 
face of the shell (classic orthotropic^^)or is so close that the eccen- 
tricity can be neglected. 

The polar orthotropic test shells, which are shells with constant 
rigidity along the main rigidity axes forming a polar coordinate system, 
have been fabricated by cementing 72 conical strips uniformly distri- 
buted over the circumference along the generatrix to an isotropic skin 
(Mylar). These orthotropic shells had the same geometrical dimensions 
and the same cross-sectional area (same weight) as the isotropic con- 
ical shell J30. 

The calculation of the rigidities of these test shells has been done 
under the following assiomptions: 

1) The reinforcements are free from the effect of transverse 
contraction. 

2) Shearing is exclusively transmitted by the isotropic skin. 

The contribution of the reinforcement to the torsional rigidity of the 

5 e 

shell is calculated according to Flugge . 

3) The flexural rigidity Du is determined for one beam element- - 
consisting of skin and reinforcement- -around the common neutral plane 
(Figure 13a). 

4) The effect of the eccentricity of the neutral plane from the 
central surface of the shell on the law of elasticity is not taken into 
account so that the elasticity equations may be written in the decoupled 
form (2. 1). 

5) It is assumed that when reaching the critical load skin and 
reinforcement buckle simultaneously. 

For the orthotropic test shell- -designation 0 30-1 --the following 
shell parameters have been determined. 

Geometry parameters: 

Zo =Zn2 = 0. 693147 \ = 2. 24471 • 10^ . 
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Rigidity characteristics: 


Cji Ci2 

7 : — = 0.519 , - 7 ; = - 0 . 154 , 

^22 ^22 

D22 Di2 

— — = 0.0573 0. 0172 , 

-L'li h)ii 


2Ci2 + C33 
2||Di2 + 2033^ 


2.41 


=0.435 . 


The results of the stability calculation together with the measured 
values are depicted in Figure 12. The straight line (I) represents the 
buckling values according to the linear theory, according to formula 
(3. 18), and illustrates properly the qualitative relationship between 
buckling load and inner pressure. The buckling loads themselves, 
however, are above the experimental values. 

The critical loads for a pin-jointed (III) and a fully restrained cone 
(IV) determined according to the nonlinear shell theory- -formulas 
(2. 21) and (2. 38)--agree with the measured buckling values for small 
inner pressure. On the one side this can be explained in such a way 
that in the reinforced test shells local prebuckling can occur only 
between the reinforcements which have respectively less effect on the 
buckling load. On the other side, the strong decrease of the theoretical 
buckling values for p = 0 is to be traced to the orthotropic rigidity 
parameters which have been calculated only approximately (especially 
the assumption that the reinforcements are free from transverse con- 
traction and that shearing is transmitted only by the isotropic skin 
represent a course approximation to the here investigated reinforced 
shells) . 

For high inner pressure, the measured values increase more than 
the theoretically determined buckling loads (curve IV for the boundary 
condition: fully restrained). This can be explained in such a way that 
due to the inner pressure the isotropic skin is stretched between the 
reinforcements (Figure 13b), whereby the axial flexural rigidity is 
additionally increased (corrugated sheet metal). The increase of the 
flexural rigidity Du with increasing inner pressure has been attempted 
to comprehend theoretically by calculating the rigidity parameter Du 
around the common neutral plane of an infinitely long segment of a 
cylinder shell with sectionally constant rigidity precurved under inner 
pressure (Figure 13b: Section I - isotropic skin. Section II - skin and 

reinforcement). The buckling loads (curve V) determined according 
to the nonlinear theory of shells taking into account the flexural rigidity 
dependence on the inner pressure agree well with the experimental 
values even for high inner pressure. 
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Figure 13. The Deformation 
Due to Inner Pressure 


CENTRAL SURFACE OF THE 



of the Reinforced Conical Shell 


A comparison between the buckling loads (Figures 10 and 12) found 
for the orthotropic (longitudinally reinforced) conical shell 030-1 and 
that for the isotropic comparison shell 130 (same weight) shows that 
under inner pressure the orthotropic shell buckles at higher loads than 
the isotropic shell. 

4.4 The Effect of Axially Symmetric Prebuckling on the Stability 
Behavior 

In the system of differential equations (2. 21) --for the deter- 
mination of the basic state of stresses--and (2. 38)--to determine the 
buckling load- -axially symmetric prebuckling Qa(^) been taken 
into account. For the approximation method (difference method) 
described here it does not present any difficulties to investigate arbi- 
trarily axially symmetrically precurved conical shells. 

Since the predeformation of the test shells is not known, a fictitious 
axially symmetric prebuckling of the following form is assumed for 
the theoretical investigations: 
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0 


for 


z - ZA 


WA 




cos 


3 


(z - ZA) 
bA 


TT 


2 


for 


z - zA 



(4.6) 


Here designates the prebuckling amplitude related to the wall 
thickness, represents the width of the prebuckling, and gives 
the position of the prebuckling maximum along the generatrix. 

The numerical evaluation of the systems of differential equations 
(2. 21) and (2. 38) results in the critical loads plotted versus the pre- 
buckling amplitude in Figure 14 for a prebuckling according to Equa- 
tion (4. 6) with bA/zo = 0* 25, za/^o = 0- 5. As was to be expected, the 
buckling load decreases with increasing amplitude ^o* 

Figure 15 illustrates the effect of the position of prebuckling on 
the critical load. It is seen that prebuckling at the edge z = 0 results 
in a greater decrease of the load than the same prebuckling at the edge 
z = zq. This tendency is even increased with inner pressure as is seen 
from a comparison of the curves II {prebuckling at the edge z = 0) and 
III (prebuckling at the edge z = zq) in Figure 16. 

An additional phenomenon is to be noted in Figure 16. 

With increasing inner pressure in agreement with experiment and 
theory, the effect of prebuckling on the critical load becomes smaller 
and smaller. The buckling values measured for four different test 
shells scatter greatly for low inner pres sure- -due to the predeformation 
which is different for each test shell--while for high inner pressure the 
experimental buckling loads are about the same and agree well with the 
theoretical buckling loads of the ideal as well as of the precurved fully 
restrained conical shell. 

Taking into account prebuckling, the decrease of the buckling loads 
can be explained for low inner pressure. An absolute comparison 
between measured and calculated buckling loads is, however, not per- 
mitted since no measured results are available concerning the pre- 
deformation of test shells and prebuckling on which the calculation has 
been based has been assumed arbitrarily. 
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Figure 14. Effect of the Prebuckling Amplitude on the 
Buckling Load of an Isotropic Cone (J30, p = 0) 




Configuration of the Predeformation: 


I 




Figure 15. Effect of the Position of Prebuckling on the 
Buckling Load of an Isotropic Cone (J30, p = 0) 







Theoretical buckling load 
according to the nonlinear 
theory (fully restraint). 



Figure 16. Comparison of the Theoretical Critical Loads of Ideal and Predeformed 
Isotropic Conical Shells (J30) Dependent on the Inner Pressure With the 
Experimental Values 






5 . Summary 


The previously known investigations of the stability behavior of 
thin-wall conical shells under axial load and inner pressure have shown 
that by means of the linear theory the qualitative effect of the geomet- 
rical characteristics and the rigidity parameters on the buckling load 
can be determined. However, the calculated critical loads deviate 
greatly from the buckling values found experimentally. To gain a better 
agreement between experiment and theory, the derivation of the dif- 
ferential equations of the state of stresses and deformation in this 
paper takes into account nonlinear members. The derived system of 
nonlinear differential equations generally describes the stress and 
deformation states of a conical shell deformed, prebuckled orthotropic 
(special case "isotropic’*) conical shell under inner pressure and 
arbitrary edge loads. 

By restricting ourselves to axially symmetric edge loads and to 
an axially symmetric predeformation, special forms of this system of 
differential equations for the basic state of stress and for the additional 
stresses and deformations at the branching point can be stated. The 
investigations in this paper are limited to these special forms of system 
of differential equations. 

The system describing the axially symmetric basic state of stresses 
consists of two coupled nonlinear ordinary differential equations of 
second order in z for the section force No and for the deformation 
parameter Qo (1st derivative of the displacement w perpendicular with 
respect to the central surface of the shell with respect to z). In the 
linear theory of stability it has been assumed that the conical shell 
assumes a state of flexure-free diaphragm stress until buckling occurs. 
It can be shown that by means of the diaphragm stress the equilibrivim 
conditions at the shell element, but not the compatibility condition and 
all boundary conditions, can be satisfied. The state of diaphragm 
stress is thus no solution of the shell problem considered here. 

A general solution in closed form of the system of nonlinear dif- 
ferential equations for the basic state of stress is not possible with 
the familiar methods. For the linearized system describing the 
stresses and deformations for small loads, a solution in closed form 
has been found for the isotropic ideal conical shell. (The particular 
solution of the linearized system is identical with the state of dia- 
phragm stress.) The iterative solution of the nonlinear system is 
performed with a difference method enabling the satisfaction of all 
boundary conditions. The finite equations at each intermediate point 
can be combined to a nonlinear matrix equation which was iteratively 
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solved by means of Newton* s method. By calculating suitable initial 
values (extrapolation with the approximate solution of next lower load 
steps) the approximate solutions according to, Newton* s method converge 
rapidly (after 3 to 4 iteration steps the agreement includes the fifth 
decimal). The convergence of the difference method is assured if the 
approximate solutions with increasing number of intermediate points 
tend toward a common limiting value. The number of the equidistant 
intermediate points was increased until the difference between two 
successive approximations remained below a specified error limit. 

The required number of points to maintain this limit increased with 
increasing load. Because of the size of the available computer (a 
maximvim of 45 intermediate points could be used in the calculation) 
the axially symmetric basic state of the stresses could not be deter- 
mined for arbitrarily large loads. For the range of inner pressures 
(0 ^p ^ 1. 0) investigated here, this load limit is above the critical 
load, i. e, , in the unstable range of the basic state of stress. 

Because of the investigation of the basic state of stress according 
to the nonlinear theory the appearance of ring buckling below the 
critical load- -as has been observed in the experiment for high inner 
pressure- -can be theoretically explained as well. 

A second special form of the general system of differential equa- 
tions can be stated for the additional stresses and deformations in the 
branching point. When reaching the critical load (branching point) two 
states of equilibrium exist which are infinitesimally close: one axially 
symmetric--identical with the basic state of stress--and one nonaxially 
symmetric characterized by the existence of nonaxially symmetric 
stress and deformation contributions in addition to the basic state of 
stresses. For the nonaxially symmetric additional stresses and defor- 
mations, a system of two linear homogeneous partial differential equa- 
tions (of fourth order for the stress function F and the deformation W) 
with homogeneous boundary conditions can be derived. The solution of 
this linear system of equations (eigen- value problem, the critical load 
represents the eigen-value) requires the -knowledge of the stresses 
and deformations of the basic state of stresses, i. e. , the solution of 
the nonlinear system for No and Qq- If, instead, the diaphragm stresses 
are inserted for the basic state of stresses, the familiar equations of 
the linear theory of stability are obtained. 

Since a general solution of this eigen-value problem in closed 
form is not known, an approximate solution has been determined by 
means of the difference method. Solving iteratively the eigen-value 
problem a smaller number of equidistant intermediate points resulted 
in the same degree of approximations as for the determination of the 
axially symmetric basic stresses and deformations. 
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The numerical evaluation of the buckling equations of the nonlinear 
theory- -nonlinear because the axially symmetric stresses and defor- 
mations determined from the nonlinear basic equations have been 
inserted in the eigen-value problem-- shows that for high inner pres- 
sure the determined buckling loads agree well with the experimental 
values. The deviation of the critical loads calculated for ideal conical 
shells and the experimental buckling values for low inner pressure 
(especially for p = 0) is to be traced to prebuckling of the test shells. 
Investigations where the effect of such predeformations on the buckling 
behavior has been taken into account are so far known only for the 
cylinder shell where mostly an affinitive prebuckling with respect to 
the buckling configuration occurring has been assumed in the calcu- 
lation, In the present paper the differential equation of the predeformed 
conical shell has been derived, and it has been solved approximately 
for the case of an arbitrary axially symmetric fixed (nonaffinitive) pre- 
buckling. 

Taking into account a predeformation, the decrease of the buckling 
loads for low inner pressure could be proved. An absolute comparison 
between the measured results and the buckling loads calculated for a 
predeformed shell is not permitted ‘since the arbitrarily assumed pre- 
deformation is certainly not in agreement with the predeformation of 
the te st shell. 

Because of the large extent of the calculation (operations with 
large matrices) the numerical evaluation was limited to a few conical 
shells. The approximation method stated can be combined to a general 
computer program with which- -assuming a suitable computer--an 
arbitrary degree of approximation can be achieved. 
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